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SUMMARY 


A theoretical analysis was made to determine how effectively a seed materi- 
al can heat hydrogen propellant in a coaxial-flow gas -core reactor. It was as- 
sumed that the seed material was spherical and that the heat-transfer mechanism 
was radiation from the hot gas core to the spherical particles and conduction 
from the particles to the surrounding hydrogen gas . 


The following ranges of conditions were used: particle diameters from 0.05 

to 1.00 micron; power outputs of the gas core from 3195 to 10 097 megawatts per 
foot of core length, which correspond to surface blackbody temperatures of 
30 000° and 40 000° R, respectively; constant density and specific heat corre- 
sponding to 500 atmospheres and 5500° R for the hydrogen propellant; and three 
values of hydrogen gas thermal conductivity corresponding to temperatures of 
2900°, 5500 , and 8000° R. 


It is concluded that the difference between the particle and bulk gas tem- 
peratures can be kept under 500° R, and thus the propellant-seed mixture can be 
treated as a continuous gray gas. Several cases of high temperature differences 
of 1500° to 2000° R between the gas and the particles are computed, but these 
large temperature differences decrease quickly away from the heat flux source. 


INTRODUCTION 

In the coaxial-flow gas -core nuclear rocket concept (refs. 1 to 3) severe 
thermal radiation problems occur. All nuclear heat released from the core trav- 
els outward toward the reactor wall. The hydrogen propellant must absorb this 
heat in order to protect the rocket wall and reach the performance levels de- 
sired in the system; however, since hydrogen gas under 10 000° R is theoretically 
transparent to radiation, the difficulty of absorbing this radiation arises. 

One of the approaches used to solve this problem is the addition of small solid 
particles to the hydrogen. The seeded hydrogen will then absorb the radiated 
heat from the core. 

In this system it is important to know how effectively the seeding parti- 
cles can transfer heat from the hot gas core to the hydrogen propellant. Since 
the primary purpose is to heat the hydrogen gas, the difference between the par- 



tide temperature and the bulk gas temperature becomes an important quantity • 

Some previous work in the area of radiant heating of gas containing parti- 
cles has been reported. Theoretical investigations of the transmission proper- 
ties of particles are given in references 4 and 5, whereas experimental data are 
reported in references 6 to 9, An analysis of temperature differences between 
a particle and a gas is reported in reference 10 for a particle in an infinite 
atmosphere. References 11 and 12 have extended this work to include a finite 
volume of gas associated with each particle. 

It is the purpose of this report to extend the present knowledge to ana- 
lyzing the radiant heating of a dispersion of spherical particles in a coaxial - 
flow reactor. In this analysis, particular interest is focused not only on the 
parameters that influence the temperature difference between the particles and 
the gas, but also on the way this temperature difference will vary with position 
in the gas -core nuclear rocket. Questions such as how does particle concentra- 
tion affect the temperature difference and what part do heat flux and particle 
size play in the analysis must be answered. If this temperature difference is 
found to be insignificant, the hydrogen-propellant - particle mixture can be 
treated as a continuous gray gas of uniform absorptivity. If this temperature 
difference is significant, a correction factor must be used in calculating tem- 
perature profiles within the cavity of the nuclear rocket. 


SYMBOLS 

A constant of integration 

a radius of particle, ft 

B constant of integration 

c.p specific heat, Btu/(lb)(°R) 

k thermal conductivity, Btu/(hr) (ft) (°R) 

L dimensionless constant, l / a 

Jo Laplace operator 

1 one-half mean distance between particles, ft 

N particle concentration, particles/cu ft 

Q radiant heat flux, Btu/(hr)(sq ft) 

q heat, Btu 

R radial position in reactor, ft 

r radial position from particle center, ft 
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S Laplacian of t 

T temperature, °R 
t dimensionless time, 

o r o 

U dimensionless variable, X(T - Tq)/Tq 

u Laplacian of U 

V velocity, ft/sec 
v Laplacian of (T - Tq)/Tq 

X dimensionless variable , r/a 

x path length, ft 
Z axial distance in reactor, ft 

a ratio of effective to actual cross-sectional area 
9 time, hr 

p density, lb/cu ft 

o' effective cross-sectional area, sq ft 

cp dimensionless constant, (pc ) /3(pc ) 

p V p g 

Subscripts : 
cl centerline 

g gas 

p particle 

0 initial 


ASSUMPTIONS 

With the gas -core nuclear reactor (fig. l) as a basis for this study, a 
model is first postulated that most clearly fits the existing situation 
(fig. 2). Various assumptions were used in forming the model so the associated 
differential equation could be solved. The assumptions were 

(l) The seeding particles are to be spherical, and their diameter is not a 
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Figure 1. - Gas-core reactor model. 



function of time. 

(2) The temperature distribution within 
a particle is neglected. 

(3) The particle absorption cross sec- 
tion is constant with temperature. 

(4) The particles receive heat only by 
radiation, reject a percentage of this heat 
to the surrounding gas only by conduction, 
and retain the rest for increasing the inter- 
nal energy of the particles. 

( 5 ) The energy radiated from the core 
surface travels outward in the radial direc- 
tion only; reradiation from the particles, 
free convection, and axial and radial con- 
duction are also assumed to be negligible. 


(6) No relative velocity exists between 

the particles and the surrounding gas. 


(7) The particles disappear (or are assumed to vaporize instantaneously) 

at a specified temperature of 8000° R, and the change in enthalpy for this proc- 
ess is zero. 

(8) No temperature jump exists at the surface of the particles since the 
mean free path of hydrogen is small compared with the diameter of the particles. 
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Some explanation should be given for several of the assumptions used in 
this study. Since a spherical coordinate system was chosen for the particle- 
gas system (fig. 2), the use of spherical particles was made to facilitate the 
ma t hemat i c s i nvolve d . 


In order to use the assumption of no temperature gradients within the par- 
ticle, the actual temperature distribution in the particle is solved for uniform 
internal absorption. For uniform internal absorption of a particle an energy 
balance on a spherical shell yields 



5aQ 

4a 


/ \ 

( pc p>p §0 


r ^ a 


CD 


where a is the ratio of radiation the particle absorbs to radiation the par- 
ticle would absorb if it were a blackbody. The maximum temperature difference 
occurs at the steady-state condition: 


k 


- #- 

Z 5r 



r ^ a 


A double integration yields 


T 


qQr 2 

8ak p 


r = 0 


r = a 


which gives 

m _ rp = 

r=0 r=a 8kp 


( 2 ) 


( 3 ) 


( 4 ) 


In the case of surface absorption the surface temperature of the particle 
is higher than that at the center. Therefore, a temperature gradient exists in 
the opposite direction from that in the case of uniform internal absorption. 

This temperature gradient approaches zero as time approaches infinity. The tem- 
perature gradient in a particle of non-uniform absorption can be assumed to be of 
the same order of magnitude as equation (4) . The maximum temperature difference 
existing in a carbon particle will range from under 1 percent to 10 percent of 
the temperature difference that exists in the gas layer surrounding the parti- 
cle. The low and high percentage levels correspond to high and low gas thermal 
conductivities, respectively. 

From preliminary work on tungsten particles of 0.1-micron diameter (ref. 

4), it appears that the absorption cross section of a particle is reasonably 
constant with temperature. With the low temperature differences found between 
the particles and the gas, relative to the particle temperature, free convection 
is neglected. Reradiation is also neglected since the temperature of the parti- 
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cle is at all times less than 30 percent of the core temperature. 


It is also assumed that no relative velocity or eddy transfer exists be- 
tween the particles and the surrounding gas. Data in reference 13, page 266, 
indicate this is true for a particle Reynolds number less than 2. From the di- 
ameter of the particle and the Stokes velocity in a gravity field, the particle 
Reynolds number for a 1-micron particle is essentially zero. In a true flow 
case where other forces besides gravity exist, such as turbulent forces, a rel- 
ative velocity larger than the Stokes velocity may exist. The increased heat 
transfer from the particle due to this relative velocity would further reduce 
the temperature difference between the gas and the particles; therefore, the 
assumption of no relative velocity creates more severe temperature differences 
between the gas and the particles than would occur in a true flow case. End 
effects are ignored so that the energy can be assumed to radiate outward in the 
radial direction only. 

The effects of particle vaporization are not included in this analysis; 
therefore, zero heat of vaporization was used to simplify the calculations. 
However, calculations were made in order to determine the maximum time for va- 
porization of the particles in question. This included the time that was re- 
quired for the particle to attain its heat of vaporization plus the heat needed 
to raise the surrounding gas layer to the vaporization temperature of the par- 
ticles. The times involved are 0.37X10"“^ to 0.37X10’"^ second for 1.0- and 0.1- 
micron-diameter particles, respectively, at a gas thermal conductivity of 
0.852 Btu per foot per hour per °R. The seeding rate was 0.314 pound per sec- 
ond, which leads to a ratio of gas sphere radius to particle radius of 94.8. 
These amounts of time are small compared with the residence time of 0.19 second 
for the particles in the reactor model itself. 

When the mean free path of hydrogen approaches the size of the particles, 
the surface temperature of the particle is higher than the immediate gas layer 
surrounding the particle. At 500 atmospheres and 5000° R, the mean free path 
of hydrogen is 0.00273 micron, which is considerably less than any particle di- 
ameter used in the analysis; therefore, no temperature jump is experienced in 
this system. 

Constant density and specific heat corresponding to 500 atmospheres and 
5500° R were used for the hydrogen propellant. Three values of hydrogen gas 
thermal conductivity were used, which corresponded to temperatures of 2900°, 
5500°, and 8000° R. The values of the hydrogen properties used were furnished 
by the author of reference 14. 


ANALYSIS 

An overall picture of the analysis will be presented before any single 
point is discussed in detail. When a sample calculation is examined, the first 
operation is computing the heat flux at the point in question. With this com- 
puted heat flux, the particle temperature is calculated from the solution of the 
differential equation controlling the model (fig. 2). At this point the total 
heat in the system and the heat contained in the particle are known as a func- 
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Figure 3. - Gas-core reactor cavity. 


tion of time. A heat balance is then used to calculate the bulk gas tempera- 
ture. 


Since the heat flux varies throughout the R,Z-plane of the reactor model 
(fig. 3), a general expression for the spatial dependence of heat flux is 
needed. In figure 3 the overall picture is shown, where at point P(R,Z) it is 
desired to find the heat flux on a particle surface. From the Beer -Lambert law 
for radiation attenuation along a linear path of length x, 

3- = exp (-crNx) (5) 

% 

where a is the absorption cross-sectional area, N is the particle concentra- 
tion, and x is the path length. As a first approximation, the absorption 
cross-sectional area a is simply the cross-sectional area of the particle 
multiplied by a dimensionless efficiency factor a; therefore a can be ex- 
pressed as 


O 

o = art a (6) 

where a is the efficiency factor, which can vary from 0 to 1 (ref. 8), and 
a is the particle radius. In this analysis, a was treated as an independent 
variable. In cylindrical coordinates, the path length x becomes &R, and in 
addition there is a geometric attenuation of heat flux that is proportional to 
the radius ratio. Therefore, for the present situation equation (l) can be 
written in the following form: 

^ exp (-aNrta 2 AR) (7) 

% R 
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This equation for Q,(R,Z) is first solved for a particular location in the re- 
actor model and is used in the equations that yield the particle and bulk gas 
temperatures . 

In order to obtain the equation for the particle temperature as a function 
of time, an analysis of the gas-particle model (fig. 2) has to be made. An en- 
ergy balance over a differential spherical shell of hydrogen gas surrounding the 
particle is performed. The resulting differential equation to solve is 


k g a 



/ ^ St 

(p Vg Se 


l > r > a 


( 8 ) 


This is the standard heat conduction equation in spherical coordinates, and it 
can be found in reference 10, page 230. 


The two boundary conditions are 


aQ(R,Z) = -4k g || + ^ (pc p ) p at r = a 

ST , , 

0 at r -l 


( 9 ) 

( 10 ) 


The first boundary condition is an equality of heat fluxes at the surface of the 
particle. The term on the left side of the equation is the heat flux received 
by the particle. The first term on the right side of the equation is the heat 
flux into the gas, and the second term is the rate of increase of internal en- 
ergy of the particle. The second boundary condition states that the tempera- 
ture gradient is zero midway between particles. The initial condition is 

T = T q at 0 = 0 (ll) 


When Laplace transformations are used to solve the differential equation, there 
is obtained the particle temperature: 

r 


T 


_ aaQ(R,Z) J 

“ J-n “ /I, \ 


4k, 


g 


exp 


3k e 


& ^ ^ ^ P 


47tahk 9 

-r - - g 

ipc p , g 


( 12 ) 


The derivation of equation (12) is given in appendix A. 


From an energy balance on the gas -particle model (fig. 2), the difference 
between the particle temperature and the bulk gas temperature can be obtained: 


T 


T =: uaQ(R, z) I 2_ _ 


4k 


exp 


g 


3k_0 


a 2 (pc p )p 


(13) 


J 
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A detailed analysis showing the steps involved in obtaining equation (l3) are 
shown in appendix B. Note at this point that the particle concentration N 
has no effect on the difference between the particle temperature and the bulk 
gas temperature. Equations (l2) and (13) were obtained by assuming L » 1, and 
for the cases used in this study L is approximately 100. 

Because of the large number of calculations, the equations were programed 
for the IBM 7091 computing machine so that a large number of cases could be 
solved in order to analyze the trends of several variables. The computer pro- 
gram used was made as general as possible for data input. The input parameters 
needed for the program are the following: 


Initial temperature of particles and gas, °R 3000 

Weight flow of particles, lb/sec 0.314 

Weight flow of gas, lb/sec 667.0 

Gas core radius, ft 1.25 

Reactor core radius, ft 5.00 

Surface temperature range of core, °R 30 000 to 40 000 

Density of particles, lb/ft^ 100 

Density of gas, lb/ft^ 0.249 

Specific heat of particles, Btu/(lb)(°R) 0.876 

Specific heat of gas, Btu/(lb)(°R) 4.95 

Thermal conductivity range of gas, Btu/( ft) (hr) ( °R) 0.387 to 2.884 

Vaporization temperature of particles, °R 8000 

Diameter range of particles, \i 0.05 to 1.00 

Ratio of effective to actual cross-sectional area 0.01 to 1.00 


From the program the output data are Tp(R,Z), T g (R,Z) , Q(R,Z), and the net heat 
contained in the gas and particles throughout the R,Z-plane. The propellant 
velocity calculated from the weight flow of hydrogen and the dimensions of the 
reactor is 36.38 feet per second. 


DISCUSSION 

In discussing the results several trends will first be shown. This will be 
done by holding all input data constant except the variable under consideration. 
This variable will assume a variety of values, and the resulting trend will be 
plotted. After these results are discussed, there will be a series of figures 
showing lines of equal T p - T g which exist in the reactor model. From these 
figures an overall picture is shown of the temperature difference between the 
particles and the gas throughout the reactor cavity. This is done for several 
* sets of different input data, and the resulting change in the equal T - T^ 

lines is shown. In addition to plotting lines of equal T p - T g , radial plots 

< of T - T are made at several locations downstream of the inlet of the re- 

P o 

actor. These radial plots illustrate the variation of T p - T g from the core 
surface out to the reactor wall. 

In figures 4 to 6 the T p - T g was obtained from the particles moving 

next to the surface of the hot gas core. Therefore, these particles receive 
the maximum amount of heat flux in the reactor model as compared to particles 
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area, a 

Figure 4. - Comparison of particle-gas tem- 
perature difference and ratio of effective to 
actual cross-sectional area. Core surface 
temperature, 40 000° R; gas thermal con- 
ductivity, 0. 852 Btu per foot per hour 
per °R. 
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Figure 5. - Comparison of particle-gas temperature differ- 
ence and particle diameter. Ratio of effective to actual 
cross-sectional area, 1. 0; gas thermal conductivity, 

0. 852 Btu per foot per hour per °R. 


moving in the reactor cavity at greater radial 
distances. Since Tp - Tg is proportional to 


the heat flux received, figures 4 to 6 rep- 
resent maximum cases for extremely high Tp - T^. 


The ratio of effective to actual cross-sectional area a is plotted 
against Tp - T g in figure 4 for several particle diameters. A significant de- 
crease in Tp - Tg is observed as the value of a decreases from 1.00 to 0.01. 

A decrease in a corresponding to a decrease in particle diameter was observed 
in reference 8 for carbon particles. From this observation if a decrease in 
T - T due to a smaller particle diameter occurs, then there will be an addi- 

Jr O 

tional lowering in Tp - T g due to the decrease in a corresponding to the 

smaller particle diameter. For figure 4 the gas thermal conductivity is 0.852 
Btu per foot per hour per °R, and the core surface temperature is 40 000° R. 

In figure 5 the variation of Tp - T g with particle diameter is shown for 

several core surface blackbody temperatures. The ratio of effective to actual 
cross-sectional area a and the gas thermal conductivity are kept constant at 
1.00 and 0.852 Btu per hour per foot per °R, respectively. 


Figure 6 shows the variation of Tp - T g with core surface blackbody tem- 
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perature for several gas thermal conductivi- 
ties. The range of gas thermal conductivi- 
ties corresponds to the maximum temperature 
range hydrogen experiences in this reactor 
concept. An increase in gas thermal conduc- 
tivity corresponds to an increase in hydrogen 
temperature; therefore, the trend is to lower 
T - T as the temperature of the hydrogen 

ir & 

increases. For figure 6 the particle diameter 
is kept at 0.50 micron, and the ratio of ef- 
fective to actual cross-sectional area a 
is kept at 1.00. 

Figure 7 presents plots of equal Tp - T g 

in the cavity of a gas -core nuclear rocket 
model for three sets of conditions. Several 
trends are shown in this figure. First, as 
the seed diameter increases from 0.2 to 
1.0 micron (figs. 7(a) and (to)), T p - T g in- 

creases markedly near the surface of the core. 
However, T - T still decreases rapidly 

ir O 

toward the cavity wall. When the core sur- 
face blackbody temperature is increased to 
40 000° R (fig. 7(c)), a continued increase 

30000 32 500 35 000 37 500 40 000 in T^ - T™ is shown. At the surface of the 

Core surface temperature, °R F & 

core T-n - T„ is very high, but again it de- 
Figure 6. - Comparison of particle-gas tempera- F & 

ture difference and core surface temperature. creases quite rapidly away from the core sur- 

Particle diameter, 0.5 micron; ratio of effec- face. Also shown in figure 7 is the rate at 
five to actual cross-sectional area, 1.0. , . n n ,, , . 

which the particles reach their vaporization 

temperature. As the particle diameter in- 
creases from 0.2 to 1.0 micron (figs. 7(a) and (b)), more than four times the 
distance or time is required to reach the vaporization temperature. However, 
when the core surface temperature is raised to 40 000° R (fig. 7(c)), the dis- 
tance required to vaporize 1.0-micron particles is just slightly longer than 
that required for a core surface temperature of 32 500° R and particles 0.2 mi- 
cron in diameter. In figure 7 the reason the lines of equal Tp - T g slope 

away from the core surface is that, after the particles start vaporizing at the 
core, the heat flux increases for any radial position as one moves in the Z- 
direction. This is due to the fact that the vaporized particles are assumed to 
be transparent. 

In figure 8, radial plots of T p - T g are given for several axial locations 
in the reactor. In this figure the rapid decrease in T p - T g is clearly seen 

as one moves outward from the core surface to the reactor cavity wall. Because 
of the extremely rapid decrease in T p - T g with increasing radial position, 

only a small percentage of hydrogen gas experiences any severe temperature lag 
behind the particles. 
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(a) Particle diameter, 0. 2 micron; core surface temperature, 32 500° R. 



(b) Particle diameter, 1.0 micron; core surface temperature, 32 300° R. 



(c) Particle diameter, 1. 0 micron; core surface temperature, 40 000° R. 


Figure 7. - Particle-gas temperature difference in gas-core reactor. Gas thermal con- 
ductivity, 0. 852 Btu per hour per foot per °R; ratio of effective to actual cross-sectional 
area, 1.00. 
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SUMMARY OF RESULTS 


The results of this study indicate some radiative heat transfer character- 
istics of a flowing gas containing spherical micron size particles. All calcu- 
lations were carried out by using bulk properties for the gas and particles , 
constant reactor dimensions, and combinations of the following physical proper- 
ties: a core surface blackbody temperature of 30 000°, 32 500°, 35 000°, or 
40 000°; a particle diameter of 0.05, 0.10, 0.20, 0.50, or 1.00 micron; a gas 
thermal conductivity of 0.387, 0.852, or 2.884 Btu per foot per hour per °R; 
and a ratio of effective to actual cross-sectional area of 1.00, 0.10, or 0.01. 
For these ranges of conditions the following trends were obtained: 

1. As the heat flux from the core increased there was a linear increase in 
the difference between particle temperature Tp and gas temperature Tg when 

all other properties were assumed constant. In the case of 0.20-micron-diameter 
particles, a gas thermal conductivity of 0.852 Btu per foot per hour per °R, and 
a ratio of effective to actual cross-sectional area of 1.00, the Tp - Tg at 

the core surface increased from 133° to 423° R for an increase in core surface 
temperature from 30 000° to 40 000° R. 

2. There was a linear increase in Tp - T g proportional to the increase in 

particle diameter. When the core surface temperature was 40 000° R, the gas 
thermal conductivity was 0.852 Btu per foot per hour per °R, and the ratio of 
effective to actual cross-sectional area was 1.00, the Tp - Tg at the core 

surface increased from 106° R for 0.05-micron particles to 2112° R for 1.00- 
micron particles. 

3. The value of T - T was inversely proportional to the gas thermal 

ir 

conductivity. For the case of 0.50-micron particles, a core surface temperature 
of 32 500° R, and a ratio of effective to actual cross-sectional area of 1.00, 
the Tp - Tg at the core surface decreased from 1013° to 136° R for a gas 

thermal conductivities of 0.387 and 2.884 Btu per foot per hour per °R, respec- 
tively. 

4. The value of T - T was proportional to the ratio of effective to 

P o 

actual cross-sectional area. For the case where the core surface temperature 
was 40 000° R, the gas thermal conductivity was 0.852 Btu per foot per hour per 
°R, and the particle diameter was 0.50 micron, the Tp - T^ at the core surface 

increased from 10° to 1056° R for ratios of effective to actual cross-sectional 
area of 0.01 and 1.00, respectively. 

5. The particle concentration was not a factor in determining Tp - Tg. 

The particle concentration, however, affected the local values of heat flux at 
the particle surface. 

6. As the distance from the heat flux source increased, the Tp - Tg de- 
creased quickly. In the case of 1.0-micron-diameter particles, a gas thermal 
conductivity of 0.852 per foot per hour per °R, and a ratio of effective to 
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actual cross-sectional area of 1.00, T p - T g decreased from 921° R at the core 

surface to 300° R at a distance of 1 foot from the core. At the wall of the re- 
actor, which is 3.75 feet from the core surface, T p - T g decreased to 31° R. 

In view of the overall effects found in this study, it can be concluded 
that the difference in temperature between the particle and the gas surrounding 
the particle is such that the propellant -particle mixture can be treated as a 
continuous gray gas. In the few extreme cases of high particle-gas temperature 
difference, a sharp decrease occurred with increasing distance from the heat 
source. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, September 28, 1965. 


15 


APPENDIX A 


DERIVATION OF PARTICLE TEMPERATURE 


In order to obtain the equation for the particle temperature as a function 
of time^ an analysis of the gas-particle model (fig. 2) has to be made. An en- 
ergy balance over a differential spherical shell of hydrogen gas surrounding 
the particle is performed. The resulting differential equation to solve is 


kgr d 

^ 5F 

With the two boundary conditions 


(> £) - 


, x St 

(p Vg 50 


Z > r > a (Al) 


5t , ^ a (p c p)p 5 t 


aQ " " 4k g 5? + 3 


5e 


at r = a 


St 

5r 


0 at r = l 


the initial condition is 


T = T a at 8-0 


(A2) 

(A3) 

(A4) 


The mean distance between particle centers l can be expressed in the follow- 
ing manner: 


i 5 = JL 

4jtN 

The following dimensionless variables are introduced: 


(A5) 


r 

X = — 

a 

(A6) 

L-i 

a 

(A7) 

y _ X (T - T 0 ) 

T 

x 0 

(A8) 

k e 0 


f— 

(A9) 

( p = p ) g a 



Rewriting equation (Al) by introducing equation (A6) yields 
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X 2 a 2 ^ 

Equation (A8) can be rewritten 


/ 2 ST\ 

( x s ) ' 


St 

Se 


t = 


UT o m 
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X o 


Taking a derivative with respect to X gives 


St 

5x 


UT o . T o Su 

X 5x 


Multiplying each side of equation (A12) by X a gives 

x2a if * - ”o a + T o " S 

The derivative of equation (A13) with respect to X is 


5x 


2 St \ 

a ttt I = - 


X a 


SU 


S 2 U 


Su 


& = - T 0 a i +T 0»^ + T oai 


or 


5 / v 2 ST \ 

5x a 5x ) 


Tq aX 


S 2 U 

Sx 2 


Introducing equation (A15) into equation (A10) results in 


X 


X 2 a 2 "° “ SX 2 


( t ° x S) ■ 


St 


The derivative of equation (A9) is 

dt = 


k. 


'S 


d0 


(pc ) a 
p'g 


The derivative of equation (A8), with X constant, is 

ST = SU — 

X 

Combining equations (A17) and (A18) into equation (A16) yields 

S 2 U SU 
SX 2 ^ 


(A10) 


(All) 


( A12 ) 


(A13) 


(A14) 


(A15) 


(A16) 


(A17) 


(A18) 


(A19) 
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(A20) 


The Laplace transform of equation (A19 ) } with the initial condition 

u(x,o) = 0 


Su(x,S) 

ax d 

The solution to equation (A2l) is 

u = A exp(- v /s’ X) + B exp(- -/s X) 


(A21) 


(A22) 


An additional change of variables is made: 


v 



= | £(U) = | [a exp (VS X ) + 


B exp 



(A23) 


The second boundary condition equation (A3) can be expressed in terms of the 
new variables as 


Xa 5X ' x 2 a - 


(A24) 


and the Laplace transformation is then 


1 dv(L,S) n 

v 2 ax 

X a 


or 

Sv(L, S) 

5x 


(A25) 


Expressing the first boundary condition equation (A2) in terms of the new 
variables by using equations (A17) and (A18) and the derivative of equations 
(A8) and ( A6 ) ^ which is 


Sr » aSX 


results in 


aQ = -4k 


t 0 au _i_ 4a (p c p) p k fi T i 


'S aX cix 


X(pc ) a 
P g 


au 
2 ai 


+ 4k 


T 0 U 

aX 2 


at X = 1 


(A26) 


(A27) 
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The Laplace transformation gives 


Let 


<xQ ^ 4k g 5v(l,S) 


ST 


ex 


^ a (pc-p)pk_ 

+ ~ — Sv (i,s) 

3 (pc ) a 

' p'g 


9 


(p c ) 

P P 

3(pc ) 
Pg 


(A28) 


(A29) 


Then 


r& + ■ ^1.=) (A30) 

g o 

When equations (A25) and (A30) are used to eliminate A and B from equation 
(A23), the result is 


v(X,S) = 


aQa 


4Xk g T 0 S 


exp^ Vs x) - f/l) expjys (2L - X)] 


|(l - Vs +9S) exp( -/S ) - ( 1 + ~t/S+ <pS) | exp 

\1 + L VS/ 


|1/S(2L - ijjl 


(A31) 


For the particle X = 1, 



v(l,S) - -28* 


4k g T Q S 1 1 + ^/g + cpS 


L 3 S 


(A33) 


For the cases used in this study L is approximately 100. Also since -v/S 
is much smaller than 1 + -yfs + cpS for all S, -^/S can he dropped from 
equation (A33) : 
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v(l,S) = + — 

4k g T 0 S \l + cpS l 3 s 

The inverse Laplace transformation of equation (A34) gives 


m m = Oga 

P 1 0 4kg 


r 






3k Q 

4jta Nk B 

- exp 


g 

, .S' 


a2 (pc p )p 

(cc p ; g 


- 


J 


(A34) 


(A35) 
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APPENDIX B 


DERIVATION OF DIFFERENCE BETWEEN PARTICLE TEMPERATURE AND BULK GAS TEMPERATURE 

From an energy balance on the gas-particle model (fig. 2), there is 
obtained the difference between the particle temperature and the bulk gas tem- 
perature. The heat contained in the gas is 


q g = rta 2 aQ0 - | Jta 3 (pc p ) p ^ 1 - exp 


3k 6 


a (P c p)p 


4naNk Q 

+ W7 


(Bl) 


which can be written as 


= $ rt(Z 3 - a 3 ) (pc V (T - T ) 


^ 3 


~P g g 0' 


(B2) 


or as 


T - T = — 
g 0 4 


JfiL 


t - a3 )(p c p) g 


(B3) 


When equation (B3) is subtracted from equation (A35) and equation (Bl) is intro- 
duced, the result is 


T - T = Jl - exp 

P g 4k s 
g 


3k e 
g 


a 2 (pc ) 


P P 


jra^qQ(9 


4jtaNk 0 


(pc P^e f t*0 3 -a 3 )( P cJ 


P'g 


+ - (pc P } P aQa 

U - a )(pc ) g 
which can be simplified to 



3k 0 

Cl- exp 

g 

a 2 (pc ) 

l 

L p p_ 


4jtaNkg0 

x ^7l 


(B4) 


T p - T g = |l + 


’(pO 


PIP 


(i 3 - a 3 )(pc p ) g 


qQa 

4k 


1 - exp 

3k g e 

4iraJMk g 0 | 

a 2 (pc ) 


l 

L P P J 

j 


3a 2 qQ0 


4 (l 3 - a 3 ) (pc ) 

\ i\v p /g 


(B5) 


Since Z 3 »> a 3 
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